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Abstract 

The primary motivation for Moyal's approach to quantum mechan- 
ics was to develop a phase space formalism for quantum phenomena 
by generalising the techniques of classical probability theory. To this 
end, Moyal introduced a quantum version of the characteristic function 
which immediately provides a probability distribution. The approach 
is sometimes perceived negatively merely as an attempt to return to 
classical notions, but the mathematics Moyal develops is simply a re- 
expression of what is at the heart of quantum mechanics, namely the 
non-commutative algebraic structure first introduced by von Neumann 
in 1931. In this paper we will establish this relation and show that the 
"distribution function" , F(P, X, t) is simply the quantum mechanical 
density matrix for a single particle. The coordinates, X and P, are 
not the coordinates of the particle but the mean co-ordinates of a cell 
structure (a 'blob') in phase space, giving an intrinsically non-local 
description of each individual particle, which becomes a point in the 
limit to order h 2 . We discuss the significance of this non-commutative 
structure on the symplectic geometry of the phase space for quantum 
processes. 

1 Introduction. 

As is well known, the Wigner-Moyal approach has its origins in an early 
attempt of Wigner [48] to find quantum corrections to the statistical prop- 
erties of thermodynamic quantum systems. A major contribution to the 
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development of the formalism that will be the subject of this paper comes 
from the seminal paper of Moyal [41] . The aim of his paper was to "refor- 
mulate the principles of quantum mechanics in purely statistical terms" . In 
doing this, he has constructed what looks like a classical (x,p) phase space 
approach to quantum phenomena even though quantum operators play a 
primary role in setting up the formalism. His aim was to find a generalisa- 
tion of the traditional methods of statistics to meet this challenge. In fact 
what he has created is a non-commutative algebraic structure based on a 
symplectic symmetry, namely, the symplectic Clifford algebra [12]. In this 
paper I will discuss how this structure arises and forms a companion to the 
orthogonal Clifford algebra that plays a central role in the description of 
the Dirac electron. Here we will discuss the role of the symplectic structure 
leaving the account of how these two algebras can be combined to a future 
paper [38] 

Moyal's work initially received considerable opposition from Dirac [40] . 
These objections and other misunderstandings have prevented the implica- 
tions of Moyal's work from being fully appreciated in the physics community 
at large although it is used successfully in quantum many body scattering 
problems [llj . In spite of this success, one is left with the belief that it is 
some form of 'semi-classical' approach, which works well in many different 
situations discussed in [11] . 

The terminology implies that it is essentially an attempt to return to a 
classical-type description. However this is not correct as it is a full quantum 
description in it own right. To appreciate this fact, it is necessary to be 
aware that the mathematical structure used by Moyal is identical to the for- 
malism developed by von Neumann in a very early paper "Die Eindeutigkeit 
der Schrodingerschen Operatoren" [42]. The content of this paper is signif- 
icantly different from that of his classic book "Mathematische Grundlagen 
der Quantenmechanik" [43] to warrant reexamination. 

In this paper I want to bring out the exact relation between Moyal's ap- 
proach [41] and this early work of von Neumann, not merely for historic rea- 
sons, but also to propose a new way forward to explore the non-commutative 
geometric aspects of the Moyal algebra that does not put the wave function 
centre stage [27] [25]. 
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2 Moyal's Proposals. 



2.1 The Characteristic Function 

To explore the possibility that the standard approach to quantum mechan- 
ics disguises a deeper statistical theory based on traditional statistical ideas, 
Moyal proposed that we start by defining one of the standard tools of a statis- 
tical theory, namely, the characteristic function. To extend these statistical 
methods to the quantum domain, he introduced a characteristic operator 
namely, the Heisenberg-Weyl operator 

M(t, 9) = e iirp+9x) (1) 

where x and p are the usual quantum operators with r and 9 being arbi- 
trary parameters. To define the characteristic function, we first take the 
expectation value of M(r, 9) by writing 

Afy(<r,0,i) = {i/>\M(t,9)\iI>) = I \p*(x,t)e^ Tp+ex ^(x,t)dx. (2) 



Using e l ( T P+ 9x ) = g iT P/ 2 e idx e lT P/ 2 ; expression (|2|) can immediately be re- 
duced to 



M^(t, 9,t)= I ij)*{x - r/2, t)e lt)x ip{x + r/2, t)dx. (3) 

Mf(r,9,t) is then defined to be the characteristic function. As in classical 
statistics the probability distribution F^(p,x,t) can be found simply by 
taking the Fourier transform of M^(t, 9,t) so that 



F^(p,x,t) = J M i ,(T,9,t)e~ i{Tp+9x) dTd9 

= ^ j iP*(x-T/2,t)e- ipT iP(x + T/2,t)dT. (4) 

This will immediately be recognised as the distribution first introduced by 
Wigner [48] . It should be noted that here p and x are continuous variables in 
some, as yet unspecified phase space and should not immediately be thought 
of as specifying the position and momentum of a single particle. Rather it 
should be thought of as some arbitrary coordinate system on a symplectic 
manifold with an invariant 2-form, a(p, x) yet to be defined^- 



1 Here we are concerned with conceptual questions so the mathematical formalism will 
be kept as simple as possible using a two-dimensional phase space. The generalisation to 
a higher dimensional phase space is straight forward. 
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Finally we can find the expectation value of any bounded operator through 
(A) = J a(p,x)F^(p,x,t)dxdp (5) 
where a(p, x) is an ordinary function on the symplectic space. 



2.2 The von Neumann Approach. 

It is interesting to note that the expression ([2]), which Moyal called a 'charac- 
teristic function', had already appeared in a classic paper by von Neumann, 
"Die Eindeutigkeit der Schrodingerschen Operatoren" [32], but von Neu- 
mann gave it no name. In his paper, von Neumann laid the foundations of 
what we now know as the Stone- von Neumann theorem, in which he showed 
that the Schrodinger representation is irreducible and unique up to a unitary 
equivalence. 

Indeed what we will now show is that the mathematical structure de- 
veloped by von Neumann is identical to the one that appeared in Moyal's 
paper. Thus the Moyal algebra lies at the very heart of the standard oper- 
ator formulation of quantum mechanics. 

Rather than starting with the well known relation 

[x,p] = ih 

von Neumann introduces a pair of bounded operators, U(a) = e iap and 
V(/3) = e^* so that the non-commutative multiplication can be written in 
the form 

U{a)V{P) = e ia Pv(P)U(a). (6) 

together with 

U{a)U{p) = U(a + /?); V{a)V(p) = V(a + 

These relations were originally introduced by Weyl |47j . 
von Neumann now defines an operator 

S(a,(3) = e- ia ^ 2 U(a)V((3) = e ia ^ 2 V((3)U(a) 

which can also be written in the form 

S(a,l3) = e ^P+m_ ( 7 ) 
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This is exactly the operator M(t, 6) introduced in equation (p}, only the 
notation is different. Thus Moyal's starting point is exactly the same as 
that of von Neumann. 

Let us go further, von Neumann then proves that the operator S(a, 0) 
can be used to define any bounded operator A on a Hilbert space through 
the relation 

A = J f a(a,P)S(a,p)dadp. (8) 

where a(r, 9, t) is the kernel of the operator. 

To proceed further, von Neumann defines the expectation value of the 
operator A as 

(il>\A\il>) = J J a(a,P){il)\S(a,0)\il))dadp. (9) 

Here 

{ip\S(a,P)\ip) = (^le^+^l^) 

so that 

(1>\S(a,P)\il>) = (iP\M(r,em, 

which, apart from a change of variables, is identical to the expression used 
by Moyal in equation ([2]). If we now use the Fourier transformation of 
M^(r,6,t) in equation Q, we find it immediately gives the Moyal equation 
([5]) for the expectation value, which for an operator that is time dependent 
gives 



(if>\A(t)\if>) = j J a(p,x,t)F^(p,x,t)dpdx. (10) 

Here a(p, x, t) is called the Weyl symbol, its Fourier transform being a(r, 6, t). 

Moyal noticed that if a(p, x, t) could be regarded as one of the possible 
values of A and if we could regard F^{p,x,t) as a probability distribution, 
then the RHS of (|10p has exactly the form of a classical expectation value. 
So why not treat F^(p,x,t) as a probability distribution? 

Unfortunately it is not difficult to find situations in which F^(p, x, t) be- 
comes negative so it cannot be a true probability. Nevertheless it can be used 
to calculate correct expectation values, hence the term "quasi-probability 
function". However the appearance of negative probabilities is not a satis- 
factory feature in a theory aimed at generalising ordinary statistical methods 
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to cover quantum phenomena. Of course Moyal was well aware of this diffi- 
culty. He was also aware that Dirac [16] had earlier discussed the appearance 
of negative probabilities in the Klein-Gordon equation, arguing that nega- 
tive probabilities should not simply be "considered as nonsense" . "They are 
well defined concepts mathematically, like a negative sum of money", an 
idea that was developed further in a paper by Bartlett [2]. 

The question of negative probabilities was later taken up again by Sudar- 
shan [45j and more recently by Feynman [18] who also argued that negative 
values were quite acceptable provided the negative values do not show up in 
any observable situation. The method can give a convenient way of calcu- 
lating expectation values which may prove more difficult to evaluate using 
standard techniques [TT] . 

In another approach Bohm and Hiley [5] have shown how these negative 
probabilities can be circumvented if desired, but this leads to a different 
theory. But notice the debate about negative probabilities only arises if we 
attempt to interpret F^{p,x,t) as a probability density in equation (fTUj) , 
Before dismissing the approach out of hand, let us look at the meaning of 
this term from two different points of view. 

2.3 Expectation values 

We have seen how the expectation value of an operator can be written in 
the form (I10p whereas the usual form of the expectation value is written as 




Let us now change to the coordinates 



x' = x - y/2; 



x" = x + y/2 



so that 




(x - y/2\A\x + y/2)p^{x - y/2, x + y/2, t)dxdy 



Now define F^(p,x,t) by 
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Substituting into equation ([lip and rearranging, we find 



<*«» = /[/ 



(x - y/2\A(t)\x + y/2)e" typ dy F^(p,x,t)dpdx. 



Defining the square bracket to be a(p,x,t), we find 



(A) 



I 



a(p, x, t)F^(p, x, t)dpdx. 



which is identical to equation ([TO]) . This result shows quite clearly that all we 
have done is to change coordinates and rearranged the standard expression 



the exact meaning of the variables (p, x). 

2.4 The Wigner Distribution and the Density Matrix. 

So far we have started from a definition of the characteristic function M^(t, 8, t) 
and introduced a pair of unidentified variables (p, x) through a Fourier trans- 
formation It is often assumed that these variables refer to the position 
and momentum of a particle but is this actually the case? 

Following on from the paper of Takabayasi [46J, Bohm and Hilcy [5 J 
showed that one can construct a symplectic phase space from a pair of 
points in a configuration space. To do this we start with a density matrix 
written in the form p(x,x'), thus regarding it as a two point function in 
configuration space. By writing this density matrix in the form 



PlP (x,x',t)=^(x,t)r(x',t) = — 0(p,t)e ixp 0*(p',t)e-™'P'dpdp' 



and introducing new co-ordinates 

X = {x + x')/2; t = x-x'; and P = {p + p')/2; 9 = p - p' (11) 
the density matrix can be transformed into 



for the expectation value of an operator A. As yet we still have not identified 




fty(X, r,t) = — f f 4>(P - 6/2, t)e i9X 0*(P + 6/2, t)e irP dPd6 




We can now write this equation in the compact form 




(12) 
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where 



F4P,X,t) = ±- J 4>*{P + Ql2,ty X6 4>(P-QI1,t)dd. 
Or taking the inverse Fourier transform we have 

F % h(P,X,t) = ^ J <P*(X -T/2,t)e- iPr iP(X + T/2,t)dr 

which, of course, is just the Wigner distribution (J3]). 

Thus we see that the coordinates used in the Moyal distribution, F 1 p(P, X, t), 
are not the position and momentum of a single particle, rather X and P are 
the mean position and momentum of what could be taken to be a cell in 
phase space. But we are using the density matrix to describe a single parti- 
cle, not an ensemble of particles. Thus our single particle is not represented 
by a point but as an extended region or 'blob' in phase space. In some earlier 
papers [31] [33] I tried to develop a quantum dynamics based on an evolution 
of a primitive cellular structure in phase space with limited success. 

To make more progress with these ideas it is necessary first to realise 
that the space we are talking about has a symplectic structure. To see this, 
let us return to equation ([7|) which will be recognised as a generator of the 
Heisenberg group. These generators satisfy 

S(a,/3)S(a',P') = e ia S(a',(3')S(a,f3) 

where a = {pi (i — a/3'), an invariant antisymmetric bilinear form which 
implies that the (a, (3) space is a symplectic space. This, in turn, implies 
that the (X, P) space is also a symplectic space which is necessary if we are 
to call the (X, P) space a 'phase space'. 

This now opens up a more fruitful approach that has been proposed by 
de Gosson [26] . It means that we are building quantum mechanics on a sym- 
plectic space, a space upon which classical mechanics is also built. Thus we 
see that classical mechanics and quantum mechanics are different structures 
built on the same underlying space, which means that this approach offers 
the possibility of giving a clearer understanding of the relation between the 
two. Symplectic spaces have topological invariants called symplectic capac- 
ities, which have shown to be intimately related to the uncertainty principle 
[22j . In this respect it seems as if classical mechanics already contains a 
shadow of the uncertainty principle. 
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2.5 The Meaning of the Weyl Smbol a(p,x,t) 

Having clarified the meaning of F^(p, x, t) H let us now examine the precise 
physical meaning of the Weyl symbol a(p,x,t) introduced in equation (|10p . 
Let us start with the usual definition of the mean value of the operator 

A, 

(A) = (V>(i)|i|^(i)) = j (iP(t)\x')(x'\A\x")(x"\iP(t))dx'dx" 

(x'\A\x")p(x' , x" , t)dx'dx" 



Let us now change coordinates using x' = x — r/2 and x" = x + r/2, then 

(A) = J (x - t/2\A\x + t/2) P (x - r/2, x + r/2, t)dxdr 
Now write p(x — r/2, x + r/2, t) = f F^(p, x, t)e ipT dp, and we find 

(A) = J {x-r/2\A\x + T/2)e' ipr dT [F i ,(p,x,t)dpdx], 
which becomes equation (jlOp if we identify 

a(p, x,t) = J(x- r/2\A(t) \x + T/2)e~ ipr dT. 

Thus we see that a(p,x, t) is a transition probability amplitude integrated 
over the 'blob' at the mean position x when the blob is moving with mean 
momentum p. The Weyl symbol a(p, x, t) is sometimes called "classical 
observable associated with the observable A", but I find that misleading 
since there is very little that is classical about a(p,x,t). It is a transition 
probability amplitude. In section (|2.7|) we will show that these symbols 
combine under a non- commutative multiplication rule, a(p,x,t) -kb(p,x,t), 
emphasising again that they are definitly not classical functions. 

Let us conclude this section by pointing out that F^(p,x,t) is the Weyl 
symbol of the density operator. To show this let us write equation (j4j) in a 
slightly different form 

F i ,(p,x,t) = ±- f e- ipT (x-T/2,t\ (IV'XV'I) \x + r/2,t)dr. 



We will follow convention and revert to lower case X and P. 
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Notice (IV'KV'I) is j us t the density operator, p, for a pure state, so that we 
can write 



Thus clearly demonstrating that the probability distribution F^(p,x,t)[:= 
p^(p, x, t)] is simply the Weyl symbol for the density matrix for a pure state 
in the (p,x,t) representation. Note there is no reason why a density matrix 
should stay positive. The positivity condition is only desirable if F^(p,x,t) 
is to be regarded probability density. 

2.6 Moyal's contribution 

Since the mathematical structure originally introduced by von Neumann is 
identical to the one later used by Moyal, I would like to clarify the signifi- 
cance of Moyal's contribution. 

von Neumann left his results expressed in terms of parameters, a, and /3, 
whose physical meanings were unspecified and so the physical significance of 
the approach was unclear. What Moyal noticed was that a simple Fourirer 
transformation will put the von Neumann expression (ip\S{a, j3)\ip) into a 
form that looked like a classical probability density in an (p, x, t) phase 
space so that equation (fT0|) became a straight forward statistical average. 
This immediately suggest that S(a, (3) could be regarded as a characteristic 
operator, a quantum generalisation of the characteristic function which plays 
a central role in the general theory of statistics. Since equation (llOj) looks 
like a classical average perhaps quantum mechanics is basically a statistical 
theory over a phase space, even though this seems to be called into question 
by the uncertainty principl^l. 

There are three things to notice about this generalisation 

• The density matrix we start with is a two point function in configura- 
tion space. 

• These two points are transformed into a one point function in the 
higher dimensional phase space, that is, a space where each point is 
now given the coordinate (p,x). 

• The density matrix is applied to the description of a single particle, 
rather than an ensemble of particles. 

3 If xi and X2 are conjugate points, then one can show [X, P] = 0. 




(13) 
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What does this all mean for the behaviour of a single particle? It suggests 
that in the quantum domain we can no longer describe the particle by a 
local entity, such as a precisely localised 'small rock' with sharp position 
x and a sharp momentum p. Its description requires a region of space in 
order to describe its behaviour. That is, it is behaving like an extended 
'blob' |26j specified by a mean position and a mean momentum. To get a 
deeper understanding of what all this means, we must explore further the 
mathematical structure of the Moyal algebra. 



2.7 The Moyal Algebra 

Let us begin by exploring the formal mathematical structure of the Moyal 
algebra a little further. What we will show is that the phase space is a non- 
commutative phase space. This has profound implications for the underlying 
symplectic geometry (See de Gosson [22], [25]). 

von Neumann showed that to completely reproduce the results of quan- 
tum mechanics, it is necessary to introduce a non-commutative star product, 
which in the (p, x) representation, is defined by 

a(p, x) * b(p, x) = J J e 2i( - pv - x ^a(p-C,x-7])b(C,r])d^drj, (14) 

where 

b(p,x)= f I e^-^m^didrj. 



Although the product was first introduced by von Neumann, it is now called 
the Moyal star product, since it was Moyal that showed it could be written 
as a formal series 



a(p, x) -k 6(p, x) = a(p, x) exp 



ih ( d d d d 
2 \dx dp dx dp 



b(p,x). (15) 



This form is extremely useful for cases where the a{p, x) and b{p, x) are finite 
polynomials. For example it is trivial to show that 

x -k p — p-k x = ih. 

showing how the famous Heisenberg defining relation appears in this phase 
space. As a passing remark, it should be noted that for non-polynomial 
functions, product rule (fT4]) should be used to avoid convergence problems. 
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The star product enables us to define two types of bracket; the Moyal 
bracket defined by 

a-kb — b-k a 

{a, o}mB = rr , 

in 

and the Baker bracket pQ (or Jordan product) defined by 

a-kb + b-k a 
{a,b\ B B = 2 ' 

Clearly the Moyal bracket replaces the quantum operator commutation re- 
lations [A, B] . 

The nature of the ^-product means that the Moyal bracket will, in gen- 
eral, be a power series in h. If we retain only the terms to 0(h), we find 

Moyal bracket — > Possion bracket. 

While in the case of the Baker bracket we find 

Baker bracket — > simple commutative product. 

Thus we find that the Moyal approach contains classical physics as a limiting 
case. There is no need to look for a correspondence between commutator 
brackets and Poisson brackets, a process which fails as is demonstrated by 
the Groenwald-van Hove theorem |30j in their well known 'no-go' theorem. 

Furthermore there is no need to introduce the notion of decoherence as 
a fundamental process to obtain the classical limit. This does not mean 
that decoherence has no role in quantum physics. It plays a vital role in 
real experiments where noise and other thermal processes enter to destroy 
quantum interference, but destroying the interference does not return us to 
the classical formalism involving Poisson brackets. 



3 Moyal Dynamics 

3.1 Moyal's Original Equation 

In order to discuss the dynamical time evolution of the quasi-probability 
density F^(p,x,t), Moyal starts from the Heisenberg equation of motion for 
M<^(t, which he writes as 

m^rAt) = r r{X)t)[ H } M]^ t)dX) 
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where H is the Hamiltonian of the system. After some working, the details 
of which can be found in [31], we find 



= t 

dt 



J [H(p + 9/2, x - r/2) - H(p - 9/2, x + r/2)] 

xF^(p,x.t)e i{rp+0x) dpdx. (16) 



By using 



M^r,9,t) = F r p(ri,Z,t)e* n **QdTidS, 



equation (JT6J) can be put in the form 



dF ^P^ t] + 2F^p, x, t) sin | fdj p - t x t p ]H(p, x) = 0. 



This equation can finally be written as 

+ {F^(p,x,t),H(p,x)} M B =0. (17) 



dt 

Thus the Moyal bracket plays a key role in the time evolution of the system. 

As we have already remarked, this equation becomes the classical Liou- 
ville equation if the Moyal bracket is expanded as a series in powers of h, 
retaining only those terms of order ti 2 . Thus the mathematical structure 
of classical mechanics emerges as an approximation when the higher order 
terms in h become negligible. However this is only a part of the story and 
we need to consider the role played by the often neglected Baker bracket in 
the dynamics. 



3.2 The Key Role played by the Wigner Density Function 

From equation we see that in the Moyal approach, it is the Wigner 
density function, not the wave function, that plays a key role. Recall this 
function is a transition probability, and not a state function, a point recog- 
nised by Groenewold [28J who called it an 'ensemble operator', using the 
word 'ensemble' even when the formalism is applied to a single quantum 
system. 

Why use a statistical concept to describe a single system? Groenewold 
was quite clear. There is a certain ambiguity in the information and control 
of a single quantum system that seems to necessitate a statistical description 
in quantum mechanics. This is a feature that has long been recognised and 
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attempts to remove this ambiguity has lead to the quest for 'hidden vari- 
ables'. Groenewold left open the question as to whether we can ultimately 
find a less ambiguous description. 

As no less ambiguous description has, in fact, been found to date, my 
own conclusion is that, in quantum processes, we have reached the limits 
of the paradigm that allows us to make a sharp, unambiguous separation 
between subject and object. In this I agree with Bohr [9j. This claim 
may come as a surprise to many since I have been deeply involved in the 
Bohm approach. It is true that this approach originally started out as 
an attempt to make such an unambiguous description but the appearance 
of an additional term, the quantum energy potential in the conservation of 
energy equation (see equation [27]), blurs this distinction as this term contains 
irreducible information about relation between subject and object, as was 
discussed in detail by Bohm and Hiley [7J. 

My own position and indeed that of David Bohm, from the sixties when 
we first became colleagues, was to realise the inadequacy of a descriptive 
form that depends on objects-in-interaction. Already in his book "Causality 
and Chance in in Modern Physics" |3j , Bohm argued that we had reached the 
limits of the mechanical paradigm and something more radical was needed. 
Indeed Bohm himself had already started to explore a process oriented de- 
scription in the late fifties, publishing a key paper on the subject in the 
early sixties [4]. I have also explained my own reasons for exploring such an 
approach elsewhere so I will not repeat them here |34j . 

Although we have motivated the Moyal approach in this paper by using 
the concept of a wave function, one can start by taking the Wigner-Moyal 
density function (j4|) as basic. But as we have already pointed out, this is a 
transition probability amplitude, the description is well suited to a process 
approach, although we will not develop this line of reasoning here. 

Again as we have shown earlier, this density function is a specific rep- 
resentation of the standard density matrix, which is a generalisation of 
the wave function, introduced to describe mixed states. Thus for exam- 
ple the more general mixed state is described by p(x, t) = c n ip n (x, t) with 
^2 \ c n\ 2 = 1- Notice, however, we have again made use of the notion of a 
wave function. But the statistical distribution function F^,{p,x,t) is poten- 
tially more general since we can start with the general defining equation 
(jlOP written in the form 



If we want to adopt this position, then we must show exactly how the wave 




(18) 
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function emerges from this more general starting point. Fortunately Baker 
PQ has already done this in a paper whose importance has been almost 
forgotten. He shows that if we use the relation 

{F(p,x,t),F(p,x,t)} BB =F{p,x,t), (19) 

then we can always write 

J e ipr F(p, x, t)dp = g*(x- r/2, t)g(x + r/2, t), 

allowing us to identify the function g(y, t) with what we traditionally call the 
'wave function'. In other words when F(p,x,t) is an idempotent satisfying 
(|19p . we can associate a wave function with the system. This function 
then enables us to use the Schrddinger algorithm to calculate the statistical 
properties of the system in the way Bohr suggested. The analogy with the 
density matrix is then quite clear since it is well known that if the density 
matrix is idempotent, it describes a pure state. However we want to stress 
that in the Moyal approach, the wave function is a derived notion, a point 
that has been stressed by de Gosson [25] . 

What is important from the point of view we are developing here is the 
idempotent nature of the density function. As we have already indicated 
above we are replacing objects- in- interaction with the notion of a structure 
process introduced by Bohm [3]. What, then, is an 'object' in a process 
philosophy? Since all is process, a system that persists like a 'particle' must 
be something that keeps transforming into itself. The flame of a candel is 
something that continues to exist, yet the gasses giving rise to the flame are 
being replaced continuously provides a nice metaphore. 

In a detailed analysis of the notion of 'existence', Eddington [17] repudi- 
ates any metaphysical concept of existence like a 'particle', and replaces it 
by a structural concept of existence which is mathematically defined by an 
idempotenlQ. We are also adopting this idea and regarding objects as idem- 
potents in the process, hence we have a straight forward way to account for 
the instability of quantum particles in general, as well as pair creation and 
annihilation. Notice also a quantum object cannot be isolated clas- 
sical object. They are a feature of the overall structure process. Without 
the process, they do not exist. Alter the overall process and the proper- 
ties of the individual object change. This then goes a long way to account 
for Bohr's insistence on a key notion of quantum phenomena, namely, the 
non-separability of observed system from the means of observation. 

4 The eigenvalues of an idempotent operator are 1 or 0, to exist or not to exist. 
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3.3 The Complete Dynamics 

Let us return to the time evolution of the density function, F^(p,x,t), but 
now assuming our system to be in a pure state being specified by an idempo- 
tent density function. We have seen how the commutator in the Heisenberg 
equation of motion is replaced by the Moyal bracket, but this does not com- 
pletely determine the the time evolution of the system as has indeed already 
been pointed out by Carruthers and Zachariasen [11] . We must find the role 
played by the Baker bracket in this evolution. 

Our algebra is non-commutative, so that we have to distinguish between 
left and right multiplication, between H(p,x) * F^(p,x,t) and F^(p,x,t) * 
H(p,x), where H(p,x) is the Hamiltonian. This means that we have two 
equations for the time development, 

H{x,p)-kF^{x,p,t) = i(2ir)- 1 J e~ iTP i)*(x - t/2)1^V>0 + r/2)dr (20) 

and 

F i> (x,p,t)*H(x,p) = i(2Tr)- 1 J e- irp ^*(x-T/2) i d t ^(x + T/2)dr. (21) 

Now if we subtract these two equations, we immediately obtain equation 
(|17p . However if we add the two equations we obtain 

{H, F} BB = i(2vr)- 1 J e- irp ty*(x - r/2)^^(x + r/2)] dr. (22) 

We have introduced a condensed notation which is well known in quantum 
field theorjU. viz, 

r^t4> = r(dt4>)-(d t 4>*)i>- (23) 

We can quickly get an idea as to what the RHS means if we choose an energy 
eigenstate, ip(x,T,t) = <p(x, r)e . We find 

i{2v)- 1 j e~ lrp ^*{x-T/2)^ t ^{x + T/2) dr = -2EF^(p,x,t). 



5 The distinction between left and right multiplication is necessary even in conventional 
quantum field theory when one deals with the Pauli and Dirac particles. The double arrow 
symbol (|23p is used for the energy term, for example, in the Lagrangian for the Dirac field 
[44| . It is therefore not surprising that equation (|22|) involves energy. 
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So clearly the Baker bracket has something to do with the energy of the 
system. Therefore let us condense the notation by writing 



i(2n) 



-i 



/ 



e~ iTp iP*{x-t/2) d t ^{x + T/2) dr 



£(p,x,t) 



so that we can write equation ()22[) in the form 



S(p,x,t) + {H,F} = 



Notice also that if we use a cross- Wigner function which employs two 
different energy eigenfunctions, then the Baker bracket measures the mean 
energy of the two eigenstates. Such an equation was first introduced by 
Dahl [14J to supplement the Liouville equation in order to obtain a complete 
specification of the energy eigenstates of molecules 

The full significance of equation ()22p is still not obvious from the form 
of the the RHS of equation ([22]) but a further insight can be found by 
going to the limit 0(fr?). After some work, including writing ip(x,t) = 
R(x, t) ex.p[iS(x, £)], we find 



From the definition of the Baker bracket we find that in this limit equation 
([22]) becomes 



This will immediately be recognised as the classical Hamilton-Jacobi equa- 
tion. Thus equation (|22p is the quantum generalisation of the classical 
Hamilton-Jacobi equation. In the next section we will relate this to the 
quantum Hamilton-Jacobi equation introduced by Bohm and Hiley [7]. 

3.4 Summary of time development equations 

Since we are dealing with a non-commutative structure, we have to distin- 
guish between left multiplication in equation (|20p and right multiplication in 
equation (j2~Tj) . These equations give rise to two time development equations, 
the first being equation (|17p was written in the form 



{H,F} BB 



2{d t S)F + 0(h 2 ) 



dF^,(p,x,t) 
dt 



+ {F^(p,x,t),H(p,x)}MB =0 



(24) 
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The second equation is 



£{p,x,t) + {H 1 F} BB = (25) 

Let us emphasise again that the Moyal bracket is equivalent to the com- 
mutator 

{H,F} M b o [#,p]_ 
while the Baker bracket is the Moyal equivalent of the anti-commutator 

{H,F} BB <* [H,f>} + 



The operator equation corresponding to equation is the quantum 
Liouville equation 



while equation (|25p has the operator equivalent is 

2d t Sp + [p,H}+ = 0. (26) 

These operator equations are obtained by taking conditional expectation 
values, i.e. integrating over p. For an alternative derivation of these operator 
equations see Brown and Hiley |10j . 

Equation (|26p is very significant for those with interest in the Bohm 
interpretation. Let us choose the Hamiltonian to be H = p 2 /2m + V, then 
equation (|26l) becomes 



8 t S(x, t) + [d x S(x, t)} 2 + Q(x, t) + V(x, t) = (27) 
2m 

where Q(x,t) = V 2 R(x,t)/2mR(x,t), the quantum potential energy. Equa- 
tion (|27p is just the quantum Hamilton-Jacobi equation, the real part of 
the Schrddinger equation used in the Bohm approach. A different way of 
showing that the equations used in the Bohm model can be derived from the 
Moyal algebra will be found in Hiley [37] . But this is not simply another way 
of deriving the Bohm results. It is vital for generalising the Bohm model to 
include relativistic spin as shown by Hiley and Callaghan [35J. 



18 



4 Conclusions. 



We have shown here that the algebraic structure of the Moyal approach [31] 
is isomorphic to the algebraic approach introduced in a much earlier paper 
by von Neumann [43j. Thus Moyal uses an algebra without modification 
that is at the very heart of the quantum formalism. However von Neumann 
offers no interpretation of the formalism at all, whereas Moyal provides a 
specific physical interpretation in terms of a statistical approach based on 
probabilities in a generalised phase space. Unfortunately the interpretation 
leads to a problem, namely, that probabilities can take on negative values 
and this has left the general impression that the Moyal algebra is inade- 
quate in some way, in spite of discussions by Dirac [16] , amongst others [2J , 
|18j . that accommodate negative probabilities. It should be noted that the 
Moyal-von Neumann algebra is sound quantum mechanics and it is only the 
interpretation that is blighted. 

We have also confirmed Baker's pQ original conclusion that the Wigner 
distribution is simply the density matrix expressed in a special representa- 
tion involving the mean position of a pair of points in configuration space. 
The Wigner-Moyal transformation then enables us to construct a six-dimensional 
non- commutative phase space for a single particle [5]. In this space we see 
that the density matrix (Wigner function) describes, not a point particle, 
but an extended structure, the coordinates (p, x) describe the mean momen- 
tum and mean position of the structure, de Gosson has called this structure, 
the quantum blob [23] , 

The 'blob' structure was hinted at by Baker pQ, who argued that the 
Wigner-Moyal formalism introduces a kind of "smeared-out" projection op- 
erator for a region in phase space. However recently de Gosson [23] [22] has 
taken this further by pointing out that this structure can be related to a rich 
topological structure that underlies symplectic spaces as demonstrated by 
Gromov's [29] "no squeezing theorem". This theorem shows that there are 
areas of phase space involving pairs of conjugate co-ordinates that cannot 
be reduced in size even under a classical symplectomorphism. This is a kind 
of classical harbinger of the quantum uncertainty principle. What quantum 
mechanics does is to introduce a minimum value, h, for this area. More 
generally these 'quantum blobs' can be discussed formally in terms of the 
notion of a symplectic capacity. Further work along these lines has been 
reported by de Gosson [23] . 

As a final remark I would like to connect the ideas outlined in this 
paper to the more radical ideas I have been pursuing elsewhere [33|,[8J. 
Strictly speaking the Wigner function is a transition probability function 
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which suggests the approach be based on the fundamental notion of process. 
Without a basic process philosophy, it would be difficult to provide a physical 
motivation for giving primary relevance to the non-commutative structure. 
It should be noted that within this structure, not only do we have quantum 
mechanics, but we also have classical mechanics arising naturally in the limit 
of order h. Thus this algebraic discription provides a natural approach to 
both quantum and classical physics so that there is no need to call on the 
notion of 'decoherence' to explain the emergence of the classical world. 
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